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I. INTRODUCTION 

Skyrme Model [l| is an effective field theory describing hadrons 

a, a 

. It is defined by 

the Lagrangian density 

C s = - 4c 2 tr [tfd lt g)tf&'g)] + |tr Qg^g^gy^g^g]) , (1.1) 

where g(x) is an element of SU{2) and c 2 and c 4 are coupling constants. If we define A^(x) 
and H* v {x) by 

^ = itr(rV^), (1-2) 

H% = d^A* - d v A% (1.3) 

£5 is expressed as 

£ 5 = 8c 2 ^^ - Ci H« u H a ^, (1.4) 
where T a (o; = 1,2,3) are Pauli matrices. By definition, A"(x) satisfies the condition 

d,AZ - d v A° = 2e a ^A?Al (1.5) 
The field equation is given as the conservation law 

= °. (1-6) 
where J a,At (a = 1, 2, 3, fj, — 0, 1, 2, 3) are defined by 

J°'^ = 2c 2 A a '" + c A e a ^H^A1. (1.7) 

The field J a,tl (x) is proportional to the isospin current of the model. Another important 
current of the Skyrme model is the baryon number current N x (x) [l] defined by 

N x = -^e^e^A^AI. (1.8) 

127T 2 ^ P 

The conservation law d\N x = follows solely from the definition of N x (x) irrespective of the 
field equation for A°(x). It was shown that solitons of the Skyrme model possess polyhedral 
structures |4|, M. As for the analytic solutions for these models, only a few simple examples 

n n □ n 

|1|, |6|, lH are known. Skyrme [H found that the configuration 

g(x) = h(k-x) (1.9) 
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with k-x = k^ leading to A*(x) = k^f a (k -x) satisfies the field equation if fc M is light-like: 

x) vanishes and g(x) is independent of the coupling 



k 2 = 0. For these solutions, the field H 
constants c 2 and c 4 . In a recent paper 
solution of the form 



Yamashita and the present authors obtained a 



g(x) = h(k ■ x, I ■ x, m ■ x), (1-10) 

with k, I, m being three momenta satisfying k 2 = I 2 = m 2 = 0. In that case, the field A®(x) 
is expressed as 

Aftx) = ^a a (u,u') + \b«{u,u') + ^fc a (u,J), (1.11) 

Kj Kj Kj 

where the variables uj and uj' are linear in x^ and k 1 is defined as 



V c 2 ■ k k ) y ' 

with the triplet (i,j,k) being (1,2,3) or (2,3,1) or (3,1,2) and k l = k, k 2 = I, k 3 = m. In 
this case, we can see that the solutions dependent on C2, C4 nontrivially, and that the field 
H™(x) and the baryon number density are nonvanishing. Under the Ansatz used in ref. 0, 
the quantities \a\, \b\ and \c\ are constants and a • b, b • c and c • a are described with the 
help of the function 

= e 3 + (63 " 6 f 63 " 62) (1.13) 
piu) - e 3 

= e 3 + (e 2 - e 3 )sn 2 (^/e l - e 3 u, _ , 

where p(z) is the Weierstrass go function satisfying the differential equation 

ip'(z)] 2 = 4 [p(z) - Cl ] [p(«) - e 2 ] [p(«) - e 3 ] . (1.14) 

Here the constants ei,e 2 and e 3 are complicated functions of \a\, \b\ and \c\. They are real 
and are assumed to satisfy e± > e 2 > e 3 . sn(w, fc) is the Jacobi elliptic function of u with 
the modulus k, 2u 3 is the second fundamental period of p(z), and uj = L ■ x is a linear 
combination of A; • x, I ■ x and m • x. We note that L 2 is equal to c 2 /c4 multiplied by a 
constant independent of the momenta k, I and m. The function X(u;) oscillates even for 
large Noting that the static energy density and the baryon number density are described 
by K(cu) and its derivative,respectively, we see that they also oscillate at the spatial infinity. 



In this paper, we consider the limit that the function sn(u, k) appearing in Eq. ()1.13|) 
tends to Tanh(u). Then the simple behavior of Tanh(u) suggests that we might obtain a 
domain wall configuration in the limit. This limit is realized if the parameter e\ tends to 
e 2 . Because of the complicated dependency of (ex, e 2 ,e 3 ) on \a\, \b\ and |c|, it is necessary 
to show that there indeed exists a set (\a\, \b\, |c|) realizing e x — > e 2 and K{u) — > e 3 + (e 2 — 
e 3 ) [Tanh(A/e 2 - e 3 ^)] 2 . 

As will be shown in later sections, this is the case. As is seen from Eq. (jl.lj) . the static 
energy density of the Skyrme model consists of two terms: the one quadratic and the other 
quartic in field variables. It turns out that the behavior of the quartic term is like that of a 
domain wall. It approaches to zero at points far apart from the wall. The behavior of the 
quadratic term is also like that of a domain wall. It approaches to a constant at points far 
apart from the wall. It should be noted that the last constant is non- vanishing The baryon 
number density concentrates near the wall. The total baryon number, however, vanishes. 

This paper is organized as follows. In|nl we briefly introduce the solutions of the Skyrme 
model described by the p function. In we explain how to explore domain wall solutions, 
the detail of the analysis is given. H^|is devoted to a summary and discussions. 

II. SOLUTIONS DESCRIBED BY p FUNCTION 
A. Field equation 

To be self-contained, we briefly review the method of ref. [f| in this section. Introducing 

fby 

C = ^ (2 = 1,2,3) (2.1) 

K 

and write £ = rj = £ 2 and ( = £ 3 , we find that the integrability condition yields 

da db 

— — =2(bx a), 

Or] d£ v ; ' 

da dc , 
g--^ =2(cxa), (2.2) 

do dc 

3C-^l =2(CX6) ' 
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Here, a,b and c are three-dimensional vectors introduced in Eq.(l.ll). If we use £,77, and 
C, the field equation (jl.6|) becomes 

where D, E, F are given by 

D = a x (b x a) + a x (c x a) + c x (b x a) + b x (c x a), 

£ = bx(oxb)+hx(cx5) + cx(ox6) + ox(cxb), (2.4) 

F = cx (a x c) + c x (6 x c) + a x (6 x c) + 6 x (axe). 

To solve this highly nonlinear field equation, we must introduce some Ansatze. 



B. Ansatze 



The Ansatze adopted in ref. |6J is summarized as 



/ da 


da 


da \ 


oh 


drt 
db 


do 


dc 


dr] 
dc 


dc 




dr] 


dC ) 



( A 7 0^ 



7 + 2 -p 




/ 



(a x b) 




o a 
a + 2 -1/ 



(b x c) 



— k (3 + 2^ 



+ 



V 




P (i 



[c x a 



(2.5) 



/ 



where a, /3, 7, k, A, /j, z/, p and a are constants. It was found in ref. |6[ that the Ansatze is 
consistent with the integrability condition ()2.2|) and the field equation (j2.3J) if the parameters 
satisfy 

p a (a + 2) 



P 



ua + (3 + 2'' 



v + (a + 2)(P + 2) 

A = r, 7 = - 



pa + /3 + 2' 

W + 2) 
p 



a + /5 + 2 
a(a + 2) 



a 



v 



(2.6) 



Therefore we can regard the four constants a,P, p and z/ as independent. 
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C. Solution of the Ansatze 



From the above Ansatze, we see that a , b and c are constants. We also obtain 



a ■ b = (a + (3 + 2)J(u) + di, 
b ■ c = pj(uj) + d%, 
c • a = —vJ(lo) + d 3 , 

\x dJ{ui) 



a ■ (b x c) = — a(a + 2) 



uj = pn£- par] + /ia(, 



v duj 



(2.7) 
(2.8) 
(2.9) 

(2.10) 

(2.11) 



where J{oj) is a function of uj and di,rf 2 and d 3 are arbitrary constants. With the help of 
the identity 



[a ■ (6 x c)]" 



a 2 (a -6) (a • c) 
(6 • a) b 2 (6 ■ c) 



(2.12) 



(c - a) (c • 6) c 2 
it is concluded that J{uS) is related to the function K{uS) in Eq. ()1.13|) by 

= -J(u) - — , 

_ 2a 2 (a + 2) 2 /x 
21 ~ ~(a + /3 + 2)z/ 3 ' 

-/i 2 a 2 - z/ 2 6 2 - (a + + 2) 2 c 2 - 2pud 1 + 2(a + (3 + 2)(pd 3 - vd 2 ) 



Z 2 



6(a + (3 + 2)pp 
As is seen from (|1.13|) and (jl,14|) . K{uj) satisfies 

~dK(u^ 2 



(2.13) 
(2.14) 

(2.15) 



duo 



A[K(u)Y-g 2 K(u)-g 3 
4 [K{u) - e x ] [K{u) - e 2 ] [K(u) - e 3 ] , 



(2.16) 



where #2, #3, ei, e 2 , and e 3 are real constants satisfying e\ > e 2 > e 3 . They are related by 

(2.17) 



1 1 

ei + e 2 + e 3 = 0, e 2 e 3 + e 3 e 1 + e x e 2 = -jff2, eie 2 e 3 = -g 3 . 



00 



du 



we find that 



If we set K(uj) = p(u + const) = p(uj + u 3 ) with u 3 = § /_ 

V (ei-u)(e2— «)(e3-u) 

is regular for real values of uj and expressed by the Jacobi elliptic function by (J1.13)) . 
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D. Solution of the field equation 



The field equation is now reduced to six algebraic equations. Three of them are solved 
by fixing dx,d 2 , d 3 as 



v 4 / 



( -C(As + Bt - Cu ) ^ 



2ABC 



A(As - Bt - Cu ) 
y -B(As - Bt + Cu ) J 

(a + (3 + 7) - 3a + (u — p) + (p — a), 
(a + 13 + 7) - 3/3 + (« - A) + {y - p), 



A 

B = 

C = (a + + 7) - 37 + (p - a) + (k - A), 
s = (/? + 7 + 2 + p- p)a 2 - vb 2 + ac 2 + 2(a + 1), 
t = pa 2 + (a + 7 + 2 + v - \)b 2 - kc 2 + 2{(3 + 1), 
Xb 2 + (a + (3 + 2 + k - a)c 2 + 2(7 + 1) 



w 



-pa 



(2.18) 

(2.19) 
(2.20) 
(2.21) 
(2.22) 
(2.23) 
(2.24) 



The other three equations are reduced to the restrictions on the parameters (a, f3, p, v). 
They can be written in the form pip 2 + 2hipu + qiv 2 + 2gi + 2fiV + r« = (i = 1, 2, 3),where 
Pi,Qi, i r ii fii 9i an d fa* are polynomials of a and (3. They were solved in ref. Two examples 
of the allowed sets are given by 



(a = p,p, v) 



J_ ^5 + 4^2 
v/2 + 2 



,1,-1 



^2.25) 



and 



(a = (3,p,v) 



5 + 4V2 



V2 



,1,-1 



(2.26) 



III. EXPLORING DOMAIN WALL SOLUTIONS 



A. Conditions to realize a domain wall 



Recalling that K(u) contains the factor (e.2 — e^) sn (^/ei — e3 u 
consider some limiting cases of this factor. From the definition 



U 



dz 



y/(l-z 2 )(l-kz 2 ) 



C2-C3 
' V ei— e3 



1 2 



we here 



(3.1) 
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of the Jacobi elliptic function, we easily find 

sn(u, 0) = sinw, sn(u, 1) = tanhu. (3.2) 

The sn — > sin limit in our case is realized by e 2 — > and leads to trivial K(u) : -f^(^) = e%. 
On the other hand, the sn — > tanh limit is realized by 

ei -> e 2 , e 3 -> -2e 2 (3.3) 

and leads to nontrivial K(u): 

K(u) = -2e 2 + 3e 2 [tanhx] 2 , x = v / 3e 2 "^. (3.4) 

From the property of the function [tanhx] 2 , we expect that a domain wall-like structure 
might appear in this limit. Since ei, e 2 and e 3 depend on the parameters a, (3, /i, v and on 
the constants a, 2 , b 2 , c 2 in a complicated manner, we have to check weather this limit can 
be realized or not. 

B. A special case 

We consider the case 

H=\, v = -l, a = P = -2-^= + y^J_^l = -0.398768 • • • (3.5) 

and assume that fc, /, m are orthogonal to each other. The conditions to be considered are 
summarized as 



(1) ei = e 2 , e 3 = -2e 2 

- g 2 = 12(e 2 ) 2 , ^ 3 = -8(e 2 ) 3 . (3.6) 

^ g\ = 27gl 

(2) e 2 > -> < 0. (3.7) 

(3) a 2 6 2 > (a • 6) 2 , b 2 c 2 > (b ■ c) 2 , c 2 a 2 > (c • a) 2 for any X - (3.8) 

We note that (a • b) , (6 ■ c), (c • a) are linear in a 2 ,b 2 ,c 2 and contain [tanhx] 2 , 92 is 
a polynomial of a 2 ,b 2 , c 2 of second order, g 3 is a polynomial of a 2 ,b 2 , c 2 of third order, 
g\ — 27g 2 is a polynomial of a 2 , b 2 , c 2 of sixth order. 
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C. Static energy density 



The energy momentum tensor is given by 

Tpv = QAa,ti^v ~ Vfiv^s- (3.9) 

Setting 

k = {k ,k , 0,0), I = (l ,0,k, 0), m = (m ,0, 0,m ), /c ,^o,^o>0, (3.10) 
we have 

K 1 = . —k , K 2 = , /— /o, = \ —TTlo (3-11) 

V C 2 V C 2 V C 2 

and 

(A , A x , A 2 , A 3 ) = -(a + 6 + c, a, b, c). (3.12) 

c 4 

The variable uj defined in Eq. ()2.11j) is more explicitly given by 

uj = L ■ x, (3.13) 
Lfi = M¥^ [-vW + 2)^ + rf(P + 2)^ - 2^(/3 + 1)^], (3.14) 

We note that u depends on the time variable x . 
In the case /i = 1, v = —1, we have 

L=(L ,L), (3.15) 

L = 2A(/3 2 + 3/3 + 1), (3.16) 

L = Awu, (3.17) 
1 



-(/3(/3 + 2),/3(/? + 2),2(/3+l)), (3.18) 



A=^4#. (3-19) 
2(/3 + l)Vc 4 ' 1 ; 

w = V2/5 2 (/5 + 2) 2 + 4(/5 + l) 2 , (3.20) 
L 2 = = 2f3(f3 + 2)(/3 2 + 6/5 + 4)A 2 . (3.21) 

Since L 2 for /3 = —2 — + V 5 ^ 4 ^ j s negative, there exists a Lorentz transformation 
(Lq, L) — > (0, L'). Since uj = L ■ x = V ■ x' does not contain the time variable a;' in the new 
system, a, b and c can now be regarded as static. For selfcontainedness, we describe the 
above Lorentz transformation in some detail in the Appendix. 



The static energy density E is now given as follows. It consists of two parts E 1 and E 2 : 
8c 2 

E=—[E 1 + E 2 ], (3.22) 
c 4 

^ =2{(v 2 a 2 + v 2 2 b 2 + v 2 3 c 2 ) 

+ (2^ 2 - l)o • b + (2v 2 v 3 - 1)6 • c + (2u 3 Ui - l)c • a| (3.23) 

E 2 = (4v 2 v 3 - 1)A 2 + (4u 3 ui - l)B 2 + (4^ 2 - 1)C 2 

+ (2 + 4uf - 4v 2 v 3 - 4v 3 vi)A • B + {2 + 4vf — 4v 3 Vi - Av x v 2 )B • C 

+ (2 + 4u|-4uiu 2 -4u2U 3 )C- A, (3.24) 

where A, B, C,vi, v 2 , and t> 3 are defined by 

A = b x c, B = c x a, C = a x b, (3.25) 
Vl = V2 = cosh 6 + ^ + ^ sinh 0, t; 3 = coshfl+ ^ sinhfl, (3.26) 

ZD ZU 

sinh = -V2 - ^ + 3 ^ + 1 (3.27) 

v/-/3(/3 + 2)(/3 2 + 6/3 + 4) 

, fl / ^ + 4^3 + 6^2 + 4^ + 2 

cosy = fWTWTW (3 - 28) 

D. Explicit solution 

To obtain a set of constants (a 2 , b 2 , c 2 ) realizing the three conditions at the top of this 
section, we adopt a simplifying assumption a 2 = b 2 = x. If we set c 2 = y, we find that g 2 
and g 3 are polynomials of x, y of second and third order, respectively. It can be seen that 
many terms in g 2 and g 3 vanish when J defined by 

J = ^ + 8{3 3 + 18(3 2 + 16(3 + 4. (3.29) 
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vanishes. We here note that (3 = — 2 — i + ^—^ is one of the solutions of J = 0. Then, 

the condition (g 2 ) 3 — 27(g?) 2 = simplifies to 

y(p x + piy) (q x 2 + 8qtxy + 8q 2 y 2 ) 2 = 0, (3.30) 

p Q = -74 - 52V2 + 22)/ 5 + iV2 + 16)/ 10 + 8^, (3.31) 

Pl = -599 - 424^ + 183)/ 5 + 4^ + 130)/ 10 + 8v^, (3.32) 

g = -175979 - 124436^ + 53907^5 + 4^ + 38118)/l0 + 8^2, (3.33) 

gi = -292919 - 207125^ + 89729^5 + 4^ + 63448 )/l0 + 8^, (3.34) 

q 2 = -272444 - 192647^ + 83457-^5 + 4^ + 59013)/l0 + 8^. (3.35) 

The only solution of Eq. (j3.30j) compatible with the conditions £>0,?/>0,g2>0 and 
#3 < is 

y = tx, (3.36) 

T = C = 0.724872 • - - , (3.37) 
where C, D and F are given by 

C = -585838 - 414250^ + 179458^5 + 4^ + 126896)/l0 + 8^, (3.38) 

D = 726708054834 + 535418248054^ - 303040027330 yjl + 4^ (3.39) 

- 214281658296)/l0 + 8v / 2 + 32820149302)/ 114 + 80^, (3.40) 

F = 4( - 272444 - 192647^ + 83457^5 + 4^ + 59013)/l0 + 8^) . (3.41) 



Utilizing the above results and together with formulae (|2.7|) . (|2.8j) and ()2.9|) . condition (|3.8|) is 
satisfied automatically. Now we turn to the static energy density E. As is seen from ()3.24|) . 
the second part E 2 of E consists only of A, B and C. Under the assumptions adopted 
here, it turns out that all of lim[ x [_» 00 A, lim[ x [_K X3 B, ]im.^^ > . 00 C are proportional to 
q x 2 + 8q±xy + 8q 2 y 2 . We then see that, under the choice y = tx, E 2 vanishes at points far 
apart from the plane x = L' ■ x = 0. This behavior is just like that of a domain wall. On 
the other hand, the first part E\ of E tends to a no n- vanishing constant as x tends to 00. 
We show the behavior of A 2 , B 2 = C 2 in Figs. 1 and 2. The behavior of the static energy 
densities E\ and E 2 are shown in Figs. 3 and 4. The behavior of total static energy density 
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E = E\ + E2 is shown in Fig. 5. We note that these figures are drawn under the assumption 
x = ^Jy in addition to ()3.36|) . yielding x = r, y = r 2 . The baryon number density in this 
case is given by 



N (x) = N[fJ /2 , 



where N and no are given by 



N = n 
n 



dK{uo) 
duo 

1 W + 2) 3 



(3.42) 

(3.43) 
(3.44) 



2?r 2 (3 + 1 

The behavior of baryon number density is depicted in Fig. 6. We find that it is positive on 
one side of the plane and negative on the other side and the integral dujN (x) vanishes. 
We note that, for instance in the x = ^fy case, we have f®^ dujN(x) = — j°° dooN(x) = 
0.0297476. 

We also note that another solution similar to the above one can be obtained also in 



the case of (2.26). In this case, we have J® oo dujN( 



x 



J °° dcuN(x) = 0.0151418 for 



x 



r,y = T 



A 2 



9 A 




jo 3 




02 

2 




/ 01 





10 



10 



X 



FIG. 1: Behavior of A 2 ( a 2 = b 2 = 0.729476 • • • , c 2 = 0.532136 • 



IV. SUMMARY 



We have explored the domain wall solution of the Skyrme model. To obtain it, we 
have considered a limiting case of the previously obtained solution of the Skyrme model 
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FIG. 2: Behavior of B 2 = C 2 ( a 2 = b 2 = 0.729476 • • • , c 2 = 0.532136 • 




x 



FIG. 3: Static energy densities E\ 

described by the Weierstrass p function. In the limit considered, the function sn(u, k) tends 
to sn(u, 1) = tanh(u) leading to a static domain wall-like solution of the Skyrme model. 
Because of the complicated structure of the model, we needed to show explicitly that the 
above limit is indeed realizable. We have shown that there indeed exists a set of parameters 
which realize the limit. The two terms constituting the static energy density of the Skyrme 
model were also investigated. The behavior of the term quartic in the field variables turns out 
just like that of a domain wall : it is non- vanishing only in the neighborhood of a plane (wall) 
in the space. On the other hand, the term quadratic in the field variables approaches to 
a non-vanishing constant at points far apart from the wall. The baryon number density 
concentrates near the plane. It is positive in one side of the plane and negative on the other 
side. The total baryon number vanishes and the topological stability of the solution is not 
maintained. To obtain solutions with non-vanishing baryon number, it would be necessary 
to obtain solutions expressed not by tanh 2 (u) but by. for example, tanh(u). 
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FIG. 4: Static energy densities Ei 
Ei + E2 




X 



FIG. 5: Total static energy 
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FIG. 6: Baryon number density 
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APPENDIX A: DETAIL OF LORENTZ TRANSFORMATION 

The energy momentum tensor ()3.9|) is rewritten as 



16c 2 2 



- (IV 



+ 



C 4 



8cl 
c 4 



k^a 2 + l^b 2 + m^c 2 + r^(a ■ b) + s^{b ■ c) 



+ u^ic ■ a) - r]^(a ■ b + b ■ c + c ■ a) 

A(k^ v B ■ C + l^C ■ A + m^A ■ B) 

+ 2r^(C 2 -C-A-B-C) 
+ 2s^(A 2 -C ■ A- A B) 
+ 2u^(B 2 -B C-A B) 

- r^(A 2 + B 2 + C 2 - 2A ■ B - 2B ■ C - 2C ■ A) 



(k 



(IU J 



5 Vtiv) 



where the following matrices have been made use of: 

(\ 1 o\ 
110 

\0 0/ 

(2 1 1 o\ 

10 10 
110 
\0 0/ 



flU J 



(1 1 o\ 




/l 





A 
















10 10 


, («v) = 











^0 0/ 




V 





0/ 


(2 1 A 




<2 


1 


A 





, (*v) = 


1 





1 


10 1 











\i 1 oy 




Vi 


1 


0/ 



(Al) 



(A2) 
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We here describe some details of the Lorentz transformation from the system with L = 
(L ,L) = (L ,-\L\u) to that with V = (0, V) = (0,eV-L 2 ). It is given by 







cosh 9 sinh 9 \ I L 
sinh# cosh# / \ — \L\ 



(A3) 



from which we obtain 



sinh^ 



-V2 



P 2 + 3P+1 



^-(3(f3 + 2)(f3 2 + 6f3 + 4) 



With the help of the matrices 



Ri 



( ±- ^(A 


72 V? 



1 



and 



R2 



^ cos0 sin (/A 

1 
^— sin0 cos0^ 



we find 



R 2 Riu = — 
w 



V2/3(a + 2) cos <p + 2(0 + 1) sin <p 



\ 



^-V2(3((3 + 2) sin + 2((3 + 1) cos 



If we set cos 6 = ^2W±2l we h ave 

T 777 ' 



R 2 Riu 



So in this coordinate system, L is given by 



/ A „\ 



VUJ 







w 



/_ ir.A 



V / 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 
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Combining the above two transformations, we find that the static energy density is given 
by Eq. (1X221) . 
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